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Abstract 

The connection between dark matter and giant black holes in the galactic nuclei 
is investigated. The joint evolution of dark and luminous matter in averaged self- 
consistent gravitational potential is considered. It is shown that the distribution of 
dark matter remains spherically symmetric even in presence of essential asymmetry 
of luminous matter in the galaxy. 

The kinetic equation describing evolution of dark matter particles distribution 
function which takes into account their dynamics in gravitational potential and 
scattering on stars is derived. It is shown that significant flux of dark matter on a 
seed black hole in the center of galaxy should appear. The growth law of seed black 
hole due to absorption of dark matter is derived. It appears that during the lifetime 
of galaxy the seed black hole should grow significantly, up to scales 10 7 — 1O 8 M . 
Brief analysis of observational data shows that the presented theory is in reasonable 
agreement with observations. 

1 Introduction 

Recently, intensive development of techniques and new methods in astrophysical obser- 
vations had lead to discovery of more than 80 giant black holes in centers of galaxies 
[1], [2]. The masses of the black holes are distributed in the range 2 • 10 6 — 3 • 1O 9 M 
(M = 2 • 10 33 g is solar mass), as is shown at Fig.l. The analysis of observational data 
revealed a connection between the giant black hole mass in the center of a galaxy and 
mass of the galaxy's bulge. 

We remind that bulge is the most dense central spherical component of a galaxy, 
containing mainly old stars. The total rotational velocity of stars in bulge is usually 
significantly smaller than velocity dispersion a, that is why the form of bulge is often 
close to spherically symmetric. The density of stars in bulge grows rapidly to the center. 
The size of bulge ~ 1 — 30kpc, it grows with the galaxy mass. The masses of giant 
black holes are approximately three orders of magnitude smaller than the masses of 
host galaxy bulges. 

According to modern knowledge, the main part of matter in the Universe is made up by 
nondissipative nonbaryon dark matter. The initial density fluctuations of the dark matter 
in the early Universe had grown up forming large-scale nonlinear gravitationally bound 
objects - galactic haloes. The important feature of these haloes is singular distribution 
of density in their centers [3] . Baryon matter trapped by gravitational field of a halo has 
settled gradually to its center and formed a galaxy. The presence of singularity in the 
center of dark matter distribution allows the compressing baryon gas to form a primordial 
black hole [4] with mass of the order of 1O 3 M . Note, that in this scenario primordial 
black hole is situated exactly in the galaxy dynamical center. Later the seed black hole 
grows intensively due to accretion both baryonic and dark matter from the bulge. The 



present paper is devoted to the analysis of evolution of dark matter and its absorption by 
the black hole in the center of a galaxy. 
The paper is organized as follows. 

In the second section, the distribution of dark matter in a galaxy and its evolution 
under the influence of the averaged self-consistent gravitational field of baryon matter 
(bulge) is considered. The problem is solved in the adiabatic approximation in assumption 
of spherical symmetry of baryon matter distribution. The role of deviations from spherical 
symmetry is analyzed. It is shown that even if the baryon matter distribution is noticeably 
asymmetric, the distribution of dark matter defined by joint action of gravitational fields 
of baryon and dark matter remains spherically symmetric. 

In the third section, growth of a seed black hole in galactic center on account of 
direct capture of dark matter particles moving in the averaged self-consistent potential is 
investigated. It is shown that the mechanism is of small efficiency in real conditions and 
cannot lead to significant growth of the seed black hole mass. 

In the forth and fifth sections, the evolution of dark matter distribution function due 
to the gravitational scattering of particles by stars is considered. The scattering lead to 
violation of adiabatic approximation. The important peculiarity of the process is that 
the mean free path length of particles is much more than the size of the bulge, hence the 
collisions are rare. Taking this into account, in the section 4 we derive the collision integral 
averaged over dark matter particles oscillations, and write down the kinetic equation 
describing the evolution of dark matter distribution function. It is shown that because 
of special features of initial distribution function, the main role is played by diffusion in 
angular momentum space. 

In the fifth section, the solution of the diffusion equation is obtained, the flux of dark 
matter on black hole is derived, and the law of the black hole's growth is ascertained. 

In conclusion, we make brief analysis of observational data showing that the growth 
of giant black holes on account of dark matter absorption should be essential. 

On the whole, the presented theory is in reasonable agreement with existing observa- 
tional data concerning giant black holes. Further development both of the observations 
and the theory and their detailed comparison are undoubtedly of great interest. 

2 Adiabatic dynamics of gravitating dark matter and 
baryonic matter. 

Observation show that over distances of the order of the horizon radius our Universe is 
homogeneous, isotropic and expands uniformly. The expansion leads to a rapid cooling of 
matter. A cold gravitating gas is unstable because of the action of universal gravitational 
forces. Growth of the Jeans instability creates regions of strong compression with dimen- 
sions much smaller than the horizon radius Rh- This is of decisive importance for the 
formation (in the CDM scenario) of the large-scale structure of matters in the Universe. 

The main role in this process is played by dark matter, which manifests itself only in 
the gravitational interaction. According to the nowadays knowledge the baryonic matter 
in the Universe contribute Qb < 0.05 and thus the gravitational processes are determined 
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mostly by nonbaryonic dark matter. It is supposed as usually that the particles of which 
dark matter is composed interact very weakly with one another and with baryonic matter. 

After recombination initial linear density fluctuations Sp are growing relatively slowly 
with time t (due to expansion of Universe) according to power law: 5p/p oc t 2 ^ 3 . As soon 
as the disturbances amplitude reaches unity Sp/p oc 1 the nonlinear processes dominate. 
Well known Zeldovich-Harrison spectrum of the initial linear fluctuations is growing with 
the wave vector k: 5p(k) oc k. Due to this nonlinear stage 5p(k n )/p oc 1 is reached first 
for the small k. Nowadays the nonlinear boundary l n oc k^T 1 oc 10 — lOOMps. Thus the 
following conditions are fulfilled for the dark matter large-scale structures /, which have 
passed the nonlinear stage of the Jeans instability: 

I „ I , u I ,,N 

r -o, — «1, - « — «1. (1) 

Here // is the free path of dark matter particles which we consider as noninteracting, u 
- characteristic velocities in the inhomogeniety. We see that the nonlinear stage of Jeans 
instability is developing in noninteracting, nonrelativistic, cold gravitating gas. 

The general system of equation, describing the nonlinear gravitational compression of 
the dark matter can be presented in a form [3]: 

Aip d = AnGp 

Here p is the density, u - velocity of the gas, and ipd - gravitational potential. The nonlinear 
process is most strongly developing near the density maximums. Let us consider the 
vicinity of the one maximum, presenting the initial density distribution p in the general 
form: 

/ x 2 y 2 z 2 \ 
p t =o(r) = Po^l-^-^-^J, u t=0 = (3) 

The solution of hydrodynamical system of equations (2), lead to singularity at the 
moment of time t — t^. After that moment in a dark matter gas the multifluid motion 
is developed. It is accompanied by density oscillation. With time t — > oo the number of 
fluids is infinitely growing. The density oscillation amplitudes ai(t) at the same time is 
going to zero: Oj(i) oc 0(l/t). As a result of this kinetic mixing of fluids the stationary 
spherically symmetric cluster of dark matter is established. The main peculiarity of dark 
matter distribution in this cluster is the existence of a singularity in its core 
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Kr- 12 ' 7 (4) 



This result was obtained in analytic theory by Gurevich and Zybin [5]. 

It should be noted that numerical calculation could not receive this result for a long 
time and showed only the plane structures of the Zeldovich pancake type. Only the 
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usage of special numerical methods allowed to make a significant step demonstrating 
the existence of the stationary spherical symmetric structures [6]. But the singularity of 
density distribution was different from (4). Further development of numerical methods 
allowed to improve this result and move on very close to the law (4) [7] . So nowadays one 
can say that the numerical calculation fully confirmed the existence of a stable spherical 
clusters with the singularity law (4). This cluster was called nondissipative gravitational 
singularity (NGS). 

The main method of analytical solution [5] was based on the solution of kinetic equa- 
tion for dark matter particles using canonical action variables / = (I R ,m), where 



I R (E,m) = y-J^ p - Mr) - ^r 2 dr, (5) 

- radial action, m is the angular momentum, E - energy, returning points (roots of the 
equation E — vpd(r) — m 2 /2r 2 = 0). According to [3] the established distribution function 
is 

f(I R ,m) = f I 1 J 8 5(m 2 -l 2 I 2 R ), (6) 

where fo and lo are constants, determined by the initial distribution (3). In particular the 
case of spherical symmetric initial density distribution (3): a = b = c, is degenerate. The 
spherical symmetry is conserved at all stage of evolution, including initial compression, 
hydrodynamical singularity, kinetic mixing. In this case the angular momentum are always 
equal to zero, consequently in (6) / = 0. In the case of weak initial asymmetry: 
(a — b) 2 + (a — c) 2 << a 2 parameter l is small, proportional to e: 

l = 0.16e, e = - J (a - b) 2 + (a - c) 2 - (a - b)(a - c) (7) 

Accordingly as follows from (6) angular momentum in this case is also small and propor- 
tional to parameter e and to radial action I R . 

We will consider now the effect of baryonic matter on the distribution of dark matter. 
Let us take into account, that the baryonic matter particles lose their energy due to 
inelastic collisions accompanied by optic emission. That is why the baryonic matter 
falls down to the bottom of the potential wells, formed by the cold dark matter. The 
luminous baryonic matter then acts as indicator of the dark matter structure identifying 
in particular the position of the centre of NGS. This produces a unique object: baryonic 
matter with a halo of dark matter. Examples of such objects are galaxies: in this case 
the evidence for the presence of a dark matter halo is provided by flat rotation curves [8] . 
Other examples of such objects are clusters of galaxies, which contain a trapped hot gas, 
again directly confirming the existence of a dark matter halo [9]. 

It should be pointed out that as the baryonic matter drops to the bottom of potential 
wells, its density rises at the central regions of the wells and can become equal to the den- 
sity of nonbaryonic matter. That will lead to the considerable distortion of the canonical 
distribution of the dark matter (4). The main feature of the problem is that the baryonic 
matter loses its energy and falls to the wells bottom slowly - at the timescales comparable 
with the lifetime of the Universe. The dark matter particles trapped in the well oscillate 
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many times. It means that the radial action Ir is integral of motion (adiabatic invariant) 
and distribution function f d is conserved and described as previously by formula (6). But 
the dark matter potential in (5) should be changed to the full potential determined 
both by baryonic and dark matter ^ = ip d + ip b . Consequently the density p d takes the 
form: 

Pd (r, t) = y^r dm 2 f dE - f ^ m) (8) 
r 2 Jo Jn E JE -^(r,t) -m 2 /2r 2 



Qe is the region, where E — \l/(r, t) — m 2 /2r 2 > 0. 

V2 



I R = — / JE-V(r,t)- m 2 /2r 2 dr. (9) 

7T Jr- 

Let us consider central part of the object where M b > M d and & ip b . Potential ip b (r) is 
determined by the mass M b {r) of a baryonic matter. M b (r) is slowly changing with time 
t due to slowly falling down into the well of a barionic component. Let us suppose, for 
example, that M b has a power-like dependence on r: 

M fe (r)ocr n , ^ocr n_1 , n > 0. (10) 

From (6), (8), (9) and (10) it follows: 

Pd = Po r (9 - 39)/16 (11) 

Thus we see that in conditions (10) the power scaling for p d exists also and just as in the 
case of purely dark matter (4), the scaling law (11) does not depend on parameters of dark 
matter initial distribution l and f . For example, if M b {r) oc r 9 / 7 both p d oc p b oc r~ 12 / 7 . 
If p b oc r~ 2 , M b oc r (isothermal sphere) and p d oc r~ 15 / 8 . We see, that dark matter 
is compressed by the gravitational action of baryonic gas. For point-like baryonic mass 
distribution could be reached maximal compression of a dark matter 

Pi oc r~ 39 / 16 . (12) 

We note that the mass distribution p oc r~ m with m > 2 must lead to the black hole cre- 
ation and the problem should be considered in more details taking this fact into account. 

We discussed above the spherically symmetric distribution of baryonic matter. In real 
conditions due to rotational motion baryonic component even in the central parts both 
of the spiral and elliptic galaxies can have a noticeable nonspherical components. Can 
they lead to a significant discrepancy of dark matter distribution in comparison with the 
case of purely spherical symmetric baryonic matter? There are some general theorems 
concerning the behaviour of systems close to integrable ones. According to these theorems, 
the average action variables lR,m, m z remain constant during particles motion in weakly 
asymmetric, slowly varying potential [10]. That is why the distribution function (6) does 
not change significantly during evolution of the slowly rotating baryonic matter [11]. 
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3 Black hole growth in the center of the dark matter 
halo (self-consistent approximation). 

The aim of the present section is the investigation of possibility of seed black holes growth 
within the bounds of self-consistent, collisionless approximation. 

The growth of central black hole leads to decreasing of dark matter particles number 
in the loss cone. On the other hand, it leads to the growth of the loss cone. As the result 
of the two competing processes the flux of a dark matter to black hole can decrease. Let 
us assume that the growth of black hole stops at total mass 

M bh = M b + M d , (13) 

where M b is an initial baryonic mass of a black hole and M d is absorbed dark mass. The 
flux of dark matter on a black hole with mass M b h is composed of particles with angular 
momentum 

m <m g = 2cr g , (14) 

where r g = is the gravitational radius of a black hole. 

In spherically symmetric potential ^ (r, t) (with arbitrary time dependence) the angu- 
lar momentum of particles does not change. That is why total mass of particles captured 
by a black hole is described by initial distribution function ff. 



N (fu M bh ) = J d 3 rdhf t (r, v)0 (m g - m) , 



where 

9(x) = 1, if x > 0, and 6{x) — 0, if x < 0. The input of a dark matter particles in 
total mass of black hole is 

M d = N(UM b + M d ) (15) 

Thus the scenario of growth of black hole mass caused by dark matter particles flux 
depends on initial distribution function /j and black holes baryonic mass M b only. If the 
equation (15) has the solution M d = M d (/,; M b ) > the mass of a black hole will not 
exceed M b + M d (/$; M b ). If the equation (15) has no solution, the seed black hole will 
grow beyond all bounds. 

As an example let us consider initial isothermal distribution of a dark matter. In this 
case the distribution function, density and potential have the form: 

mrn) = — E ^e~^ 1 Po = (16) 



P(r) = §, (17) 



^(r)=47rGp ln(r), (18) 
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where aa is velocity dispersion of dark matter particles. Density of particles with angular 



momentum smaller than m g is 
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Pg = — j dE ] dm -j=±—9\E-^ d -—y{m g -m) 



or, after calculations 



fl,(r) = £|l- e -=£|, 



Hense total mass of particles with the angular momentum smaller than m g is 

n \ 1/2 



iV, = 2tt ^) m g . (19) 
Taking into an account (19) one can rewrite the equation (15) in the form: 

Md Q(l + %), (20) 



M b ^ V M b y 

where 

Q = ^(G P0 ) 1/2 
c 

The positive solution of the equation (20) exists only at Q < 1. Thus, scenario of growth 
of a black hole in the case of isothermal distribution of a dark matter does not depend 
on black hole baryonic mass but depends strongly on the dark matter density determined 
by factor p . If p is anomaly big: p > G _1 (c/87r) 2 the initial isothermal distribution 
is unstable: any small seed black hole will result in absorption of all dark matter halo. 
At p < G _1 (c/87r) 2 distribution is steady. The absorbed dark matter mass appears 
proportional to baryonic mass and is determined by expression (20). 

For our Galaxy Q ~ 0.01 and is less that 1%, what means, that the accretion of 
dark matter on the initial black hole M b lead to a very small growth of its mass only. 

Actually, distribution function of dark matter (6) differs essentially from isothermal 
one. In this case total mass of particles with angular momentum smaller than m g is 
determined by expression 

N g = (2tt) 3 / dI R / dm / dmJoI^S (m 2 - / 2 J 2 ) 9 (m g - m) , 

JO JO J-m v 7 



hence 



The equation (15) takes the form: 

M b 

where Q = Q'M^ S , Q' = (27r) 3 |/ (j^) ^ an d fo can be calculated by expression 
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M H = (2tt) 3 jf dI R | o °° dm 2 f Q lf5 (m 2 - Z 2 / 2 ) , 

where the radial action of particles near the halo's boundary is related to the halo's full 
mass M H and radius Rh by expression I max w ^G^M^rK 2 , hense 

Q ^9/8 M 7/16 

J0 8 (27r) 3 (Gi?H) 9/16 ' 

Using this relation, one can show that in our Galaxy the value of parameter Q is 
again small (of the order of 10~ 2 ), and mass of dark matter captured by black hole is not 
significant. 

We would like to pay attention that investigation of this problem in previous papers [12] 
results to opposite statement: the initial black hole was growing without any limits. The 
reason of this discrepancy is in the fact that authors did not take into account the changing 
of dark matter particles distribution function caused by black hole absorption. That is 
why the loss cone was filled constantly and the flux on black hole was not exhausted. 
Considering the transformation of distribution function one comes to our result. 

Thus, dark matter particles cannot cause a noticeable growth of the black hole because 
of conservation of their angular momentum in the spherically symmetric potential. As it 
follows from the remark in the end of Section 2, weak asymmetry of the potential does 
not change significantly this conclusion. The simplest process that can fill the loss cone is 
collisions of dark matter particles with stars. We shall discuss this point in next sections. 



4 Kinetic equation 

Gravitational interaction of dark matter particles and their collisions with stars is de- 
scribed by distribution function /(r, v,t), which satisfies the kinetic equation [13]: 

^ + {H , f} = St[f], (24) 

where Hq is Hamilton function corresponding to the motion of particles in the averaged 
self-consistent potential, St[f] is the collision term. Gravitational interaction of dark 
matter particles with individual stars is described by Coulomb law, hence the collision 
term may be written in Landau form [14] 

Si|/I = t^tj- (25) 

Here Wk p = 2nG 2 M Q A J d 3 v'wk p F(y' , r), Wk p = {u 2 5k P — UkU p ) /u 3 , Vk are the dark mat- 
ter velocity v components, F(v', r) is the stars distribution function, u = v' — v is the 
relative velocity of a star and dark matter particle, A ~ 10 is the gravitational Coulomb 
logarithm. 

In fact, as we have already noticed, the frequency of collisions of particles with stars is 
much less than the frequency of their orbital motion. Thus the kinetic equation (24) may 
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be significantly simplified. For this purpose, let us rewrite it in the action-angle variables 

% + »*wr ShM (26) 

Taking into account that the initial distribution function (6) does not depend on angle 
variables 0, we note that its time variation is determined by the collision term (25) only. 
However, the collisions are rare and St[f] oc uf, where v is the collision frequency. Hence 
we search a solution of the equation (26) in the form f = fo + vfi, where vf x is a small 
correction to f , and both the terms f fi depend on "fast" t = t and "slow" ti = vt 
times. In zero approximation by the parameter v we obtain 

dfo , df 
dt d(f)k 

From this equation follows that the main part of distribution function f does not depend 
on angle variables and "fast" time. The first approximation gives 

/ dfo dfi dfi \ r . 

u {w 1 + dr +UJ w k ) =stlM - 

Averaging this equation over the angle variables and "fast" time, taking into account 
that the second and third terms become zero, we obtain equation 

^f^ = Si[f(I,t)] (27) 

(Here we omit the index "0" in distribution function and returned to usual time t.) 

So, under assumptions on collision frequency v made above, it is possible to find the 
distribution depending on action variables only, as a solution of kinetic equation (27) with 
time-averaged collision term: 



St[f] = j^y 3 J WStiM (28) 

The time-independent kinetic equation in the averaged form (27) describing the dis- 
tribution of stars in the vicinity of a giant black hole was first studied in [15], [16] and 
[17]. A purely Coulomb potential field was considered in these works. 

Our aim is to examine time-dependent solutions of the kinetic equation (27) describing 
dynamics of particles in arbitrary central symmetric field. First, it is necessary to obtain 
the averaged collision integral (28). To do this we will follow the technique, developed in 
[18]- 

Let us consider tensor differential form (25) formally in six-dimensional phase space 
X = (v,x), assuming coordinate components of the tensor W being zero: 

Odd d 
-W kp — = -^nr-W^-^r; k,p= 1,2,3; fi, v = 1, . . . , 6 



dv k dv k dX^ * dX u 
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In this six-dimensional space we make canonical transformation to the action-angle vari- 
ables 

i-> Y,,!\ Y = {I R , m, m z ; <p u fa, 3 } • 
The differential form then takes the form 

9 w 9 1 9 rv 9 p dY »' dY »'w 

where y 7 ^ is the Jacobian. However, Jacobian of canonical transformation is unity, hence 
the differential form in new variables preserves in the form: 

d d d d 

-VV IW —— — T— 



dX^'^dX, dY,, ^ dY v , 

We search for the solution of kinetic equation (27) which does not depends on angle 
variables 0, hence -§yR-§^, = 0. After averaging over terms -§^R-§j is vanished, and the 
expression for averaged collision integral becomes again three-dimensional: 

311/1 = ik^w/' (29) 

where 

^(^/^^ < 3 °> 

To give concrete expression to the differential (29), one should know the distribution 
function of stars F. The observations of stellar dynamics in the bulge indicate that in 
the first approximation the distribution function of stars could be assumed isotropic, 
i.e. depending on energy E' — ^- + \P(r) only. For simplicity we also assume that the 
dependence of stellar distribution function on energy is power-like: 

F(v',r)=F E'-P (31) 

One can show that the potential produced by stars distributed in accordance with (31) 
is also power function of r: 

*(r) = ^ r a , a = 4/(2/3- 1) (32) 
We note that both distribution function (31) and potential (32) can be determined if 



the velocity dispersion of stars a = y {v\) = y (v 2 /3) is known. The dispersion could be 
choosen as a function of distance to the galaxy center: 

a(r) = a^l\ (33) 

The distribution function (31) and the potential (32) are uniquely dependent functions of 
the velocity dispersion, and hence the parameters (To, a only: 

2 JZ°x 2 (l +x 2 /2Y (3 dx „ a(l + a) t2 /* 

° f™x 4 (l+x 2 /2)- f3 dx {4n) 2 GJ™x 2 {l + x 2 /2)~ l3 dx 
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Note that in the case a = the relations formally lose their mathematical mean- 
ing. However, one can show that when a — > the stellar distribution function becomes 
isothermal: 



*(r) = 2cr 2 ln(r), 

and the velocity a dispersion does not depend on distance. As a rule, the parameter a in 
a bulge is rather small. 

Calculations show that for isotropic stellar distribution function F(E') tensor Wk p in 
(25) takes the form 

W kp = A(E,r)5 kp - B(E,r)^, (34) 



where 



1 fivr 2 r°° ( 1 E < E' 

' G 2 M Q A / dE'F(E') < 3v , f v , 2) ' 

A - B = ——G 2 M Q A / dE'F(E') 

3 ° Ah v ; 



3 

16tt 2 . , , . . p . .. I I . E < L' 

/*(r) 

We notes, that the radial action I R of particles moving in the potential (32) could be 
with good accuracy approximated by the expression 

I R (E,m)&J(E)-b a m, (35) 

where b a ~ 1 is some positive constant, 

J{E) = ^=e^, a = 0. 

'IT 



(We notice that in the case of Coulomb or oscillator potentials the equality (35) becomes 
exact and the constant 6=1. In the case of isothermal stellar distribution 6 pa 0.6) 

The initial function (6) does not depend on the variable m z , hence one can search for 
solution of the equation (29) as a function of I R and m. Calculating the coefficients of 
quadratic form (30) and taking into account (34), we rewrite the collision term in a form 

St[f\ = ~—m (r22— + Ri2—) + — (Ri2— + Ru— ) (36) 
m dm \ dm dI R J dI R \ dm dI R ) ' 

where 

Ru = p - 2b a q + 6 2 s, 
R\2 = q- 6 a s, R 22 = s, 



11 



2 



,i 2 -^2 

Ar -to 

where < . . . ></,= T ^m) fr- • •)> anc ^ T = 2 J r r + * is the period of dark matters partical 

radial oscillations, t> r = V2~(E — \&(r) — ^ is radial velocity. Taking into account 
that angular momentum of dark matter particles are small, in the first approximation we 
shall calculate the coefficients R a b(lR,m) in (36) at to — 0. Further, from the initial dis- 
tribution function (6) one can see that for sufficiently large time, as far as the distribution 
function has sharp maximum at to — IqIr, Iq <C 1, one can neglect the diffusion by Ir 
compared to diffusion by to. Summarizing, we rewrite the kinetic equation (27) as the 
diffusion equation depending on Ir only parametrically: 

df(I R ,m,t) , n l d d 

— = R(I R ) — —m—f{lR, to, t) 37 

at m dm dm 

Here the diffusion coefficient R = R22\ m =o = {Ar 2 )^ calculated with account of (33), takes 
the form 

R(I R ) = OAQGMQAaf^lf^ (38) 

We note that in the case of isothermal distribution of stars in the bulge a — 0, and 
the diffusion coefficient R does not depend on Ir. We emphasize that the fact that Ir 
comes into the final equation only as parameter is conditioned by the particular form of 
initial distribution function (6) containing small parameter Z - 



5 Flow of dark matter on a black hole. 

As a result of direct capture of particles in the loss cone of seed black hole, the distribution 
function of dark matter particles takes the form 

f(I R , m, 0) = f I R /s 5 (m 2 - l 2 I 2 R ) 9 (m - m g ) . (39) 

It differs from (6) by factor (to — m g ). 

The diffusion of dark matter particles in angular momentum space, their drift into 
the loss cone (14) and further capture by black hole can be described by the diffusion 
equation (37) with initial condition (39) and boundary condition 

f\m=m g = 0. (40) 

The solution of the diffusion equation can be presented in a form 

f(I R ,m,t)= dm 1 G(I R ,m,m 1 ,t)f(I R ,m,0) (41) 

Jm g 

where G is the Green function of the boundary problem (37)- (40) 
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G= dXmxZx (mi,m,) Z x (to, to 9 ) e - Ai? ( 7 «)*, 

*/ 

Z is the orthogonal system of fundamental functions of the boundary problem (37)- (40) 

J (VXm g ) N (y/Xm) - N (V\m g ) J (^/\r 



Z\ {m, m g ) 



(j 2 (y/\m g ) + N$ (yXm, 



1/2 



Jo, N are Bessel and Neumann functions of zeroth order. 

Let us now calculate the flux of dark matter on a black hole. Let D be the area in 
the phase space where Ir > 0, to > m g , —to < to 2 < to. The whole mass of dark matter 
is N(t) = J D d 3 Id 3 4>f(lR,m,t). Now, from the condition of mass conservation it follows 
that the flux through the boundary is S — —dN/dt. Taking into account (25) we obtain 



S(t) = -f d 3 Id 3 <P—R kp —f(I R ,m,t). 
Jd oh ol„ 



Let us transform now the integral over volume into an integral over surface according 
to Stokes formula. Taking into account that the only flow of matter is through the surface 
to = m g , we find that 



S(t) = 2(2tt) 3 J o dI R R(I R )m g ^/(/*, m, t) 



(42) 



Substituting the solution of diffusion equation (41) into (42), we obtain 



7T JO JO 



9 ) 



(J 2 (V\m g ) + iV 2 (y\m g )) 



1/2 : 



where 

7 = a/ (2 + a),R a = 0A6AGM Q a^. (43) 
It is convenient to introduce a normalized time T = l^ 1 R a t and make a change of variables 

2 1 1_ 

A i-> i] = \T 2 -i , Ir^ y = IqIrT 2-t , m g i-> x = m g T 2 -i 
The expression for the flux takes then a form 



*C) - ^-^^ 7 (x), (44) 



7 — 87 f°° 
( = 16 -87 ' = JL dyH ^ x ' y ^ 

•/ 



^0 (y/VX) 


N 




-N 




Jo 




Jo (V^ X ) 
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One can show that H^(x,y) is finite positive function, it is equal to zero at y = x and 
exponentially small at y > x + 4. From this follows that at the moment t the flow of dark 
matter to the black hole is made up by particles with radial action range in an interval 

HhL <lR< HhL + 4 T i/(2- 7 )_ (45) 
to to to 

Thus, the area from where the dark matter flows into black hole, grows with time as 
*V(2-7). For function $ in expression (45), there is an estimate <3> 7 (x) « rrs + 2. Finally, 
we find the expression for the flux of dark matter onto a black hole through the boundary 
of the loss cone m = m g : 

S(t) = 2<2f^^. ,46) 

The expression (46) has been obtained in assumption that the boundary of the loss 
cone m g does not depend on time. In reality, the boundary moves in accordance with 
the whole black hole mass. However, from (46) follows that the flux depends only weakly 
on the value m g . It means, that the growing of m g is a quasistationary process and one 
can use the expression (46), assuming t as a parameter: m g (t) = AGM bh (t)/c. Hence, the 
growth of the black hole is described by 

dM bh _ (2tt£ f nl/2 M* + * 

dt 7T \ C J ,1+7 01 t 1 /* • 

Assuming the seed black hole mass small, we obtain the solution of this equation 

M bh = Ct^, (47) 

To estimate the parameter f one can use the relation (23). Thus, the growth of black 
hole at the expense of absorption of dark matter scattered by stars follows a power law 
M bh oc t a , 4/7. 

6 Conclusion. 

In conclusion, we compare roughly the results of the theory presented above with the ob- 
servational data. Fig. 2 represents the dependence of stellar velocity dispersion on distance 
to the galactic center for our Galaxy and galaxies M31, NGC4258 [19]. One can see that 
the distribution of stars in bulges of M31 and NGC4258 is close to isothermal, i.e. up to 
the area of action of black hole (r < (4 — 7)pc) the velocity dispersion is almost constant, 
it does not depend on r: a = a® ~ 200km/s. In the case of isothermal distribution of 
stars in the bulge we can rewrite the expression (47): 



14 



Mbh MMo 1 (48) 

Let us suppose that parameter l of the dark matter haloes around these galaxies is 
close to average value < l >= 0.1. Implying, that masses M H and sizes Rh of haloes 
are roughly equal and constitute 1O 12 M and lOOkpc respectively we obtain from (48) 
the mass of the black hole M bh w 8 • 1O 7 M (We assume, that the age t of black hole is 
comparable with the age of the Universe t ~ 3 • 10 17 s,). Thus, theoretical predictions of 
masses of black holes in these galaxies are close enough to the observed values: M bh = 
(2.0 - 8.5) • 1O 7 M for M31 and M bh = (3.8 - 4.0) • 1O 7 M , for NGC4258. 

The central part of our Galaxy's bulge does not satisfy the assumption of isothermal 
distribution of stars at r < lOOpc. At the same time, as it is clear from (45), the dark 
matter accumulates mainly from the central region. The dependence of stellar velocity 
dispersion on distance is shown at Fig. 3 [20]. One can see that at the scales from lOpc 

to lOOpc the dispersion can be roughly represented as a = 60km/s(y^) . From (33) 
then follows a rs 0.6. In the case, the formula (47) gives for the black hole mass a value 
2 • 1O 7 M , which is much more than observed value M b h ~ 2.6 • 1O 6 M . To obtain more 
precise estimation of the black hole mass, more information on the stellar and dark matter 
distribution functions and their evolution is needed. Also, the contribution of the black 
hole into the joint gravitational potential may be important. However, we see that even 
the rough calculation of absorption of dark matter only may give a reasonable estimate 
of black hole mass. This estimate is consistent with observed masses of most giant black 
holes (see Fig.l). The fact may be considered as an additional prove of the general theory 
of the role of CDM in the large-scale structuring of the Universe [5]. 

Note also that the kinetic theory developed in Section 5 allows to describe the absorp- 
tion of luminous (baryon) matter as well as dark matter. More precisely, it is possible to 
take into account the capture of stars by black hole because of their gravitational scatter- 
ing. This process is especially important for giant black holes having mass of the order 
M b h ~ 1O 9 M and active galactic nuclei. 

Another important process described by the kinetic equation (27) is the possibility 
of diminishing of amount of dark matter in bulge. The change of dark matter particles' 
energy as a result of gravitational scattering by stars leads to their "heating" and pushing 
them out of bulge - the phenomenon analogous to Fermi acceleration. Therefore, the 
density of dark matter compared to that of luminous matter decreases. Qualitatively this 
corresponds to the results of recent observations discussed in the literature [21], [22]. 

Further development of the consistent kinetic theory discussed above and its detailed 
comparison with observational data would be undoubtedly useful for understanding of 
main physical processes in the galactic nuclei. 

The authors are grateful to V.A. Sirota and M.I. Zelnikov for numerous and useful 
discussions. 
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